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Exploring geometrical constraint effects of Cooper pairs 

Neal E. Staley and Ying LiiQ 
Department of Physics, The Pennsylvania State University, University Park, PA 16802, U.S.A. 

(Dated: March 30, 2011) 

We report low-temperature electrical transport measurements on filled Al squares prepared by 
e-beam lithography featuring a sample size ranging from 130 nm to 530 nm. The values of the 
superconducting coherence length, <f(0), were found to range from ^^100-210 nm. We found that 
phase diagrams in the magnetic field and temperature space feature quantized steps in samples 
with a size comparable with ^(0), suggesting the emergence of quantized states in these smallest 
superconducting structures. These quantized states, which are not anticipated in Ginzburg-Landau 
theory, may be associated with the quantum size effects of Cooper pairs. 



Shortly after the publication of the Bardeen, Cooper, 
Shrieffer (BCS) theory of superconductivity, Anderson 
pointed out that superconductivity is fully suppressed 
when a superconducting grain becomes so small that it's 
electron energy level spacing approaches the bulk super- 
conducting energy gap [1 , a phenomenon known as the 
quantum size effect of superconductivity. Experimen- 
tally, the superconducting energy gap of Al grains were 
found to indeed fall rapidly when the diameter of the 
grain was reduced to a critical size of Lc (material depen- 
dent, ^6-10 nm for Al)^.. The value of Lc was found to 
also depend strongly on whether the number of electrons 
in the grain is odd or even [3 . The interesting question is 
how Cooper pairs will respond to geometrical constraints 
in a grain when its size is larger than Lc but smaller 
than or comparable with the size of a Cooper pair, which 
in BCS theory is the superconducting coherence length, 
^0- The corresponding length in Ginzburg-Landau (G-L) 
theory is ^(0), which characterizes the size of the nor- 
mal core of an Abrikosov vortex and is the same order of 
magnitude as ^o- In the presence of disorder, the Cooper 
pair size is reduced - ^o becomes (^oO^^^ ^^^ i^ related 
to ^(0) by ^(0) = 0.74(^oO^^^ where / is the mean free 
path. Geometrical constraint effects of Cooper pairs my 
be expected in these ultrasmall superconductors. 

Over the past twenty years, much work on small su- 
perconductors was focused on samples with a size sev- 
eral times of ^(0), and thus much larger than the size 
of a Cooper pair. In singly connected samples such as 
filled squares or disks, in the presence of an external 
magnetic field, their behavior is determined by few vor- 
tex physics[4-7 . More recently, superconductors with 
a size smaller than or comparable with ^(0) were stud- 
ied primarily by a scanning tunneling spectroscopy probe 
[SUTO], It was found that in ultra small singly connected 
samples no vortices can exist below a critical size. Also, 
an unexpectedly large pseudogap related to modified 
electron-phonon interaction exists, as well as a quantum 
size effect enhanced superconducting gap. In addition, 
doubly connected samples have been studied, showing 
that the fundamental fluxoid quantization, which is re- 
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sponsible for the well-known Little-Parks oscillations (TT]. 
leads to emergence of a destructive regime for ultrasmall 
superconducting loops with a diameter d less than ^(0), 
as predicted by de Gennes [12^ . Here the superconducting 
transition temperature of a loop with its circumference 
less than 7r^(0) is suppressed to zero near half-integer flux 
quanta[12 , as confirmed experimentally [13 . In all these 
cases, however, the issue of the geometrical constraint 
effects of Cooper pairs was not addressed. 

In the present work we explore possible quantum size 
effects of Cooper pairs in a singly connected supercon- 
ducting square with its size similar to ^(0). Within G- 
L theory, a phenomenological framework describing su- 
perconductivity, no vortex and, furthermore, no spatial 
variation in the superconducting energy gap is allowed in 
isolated ultrasmall superconducting samples of this size 
under any applied magnetic fields. Therefore no quan- 
tized states in the superconductor are expected in G-L 
theory. Consequently, the phase diagram in magnetic 
field vs. temperature (H-T) should be featureless [14 . On 
the other hand, if we view the singly connected sample 
as being composed of a series of doubly connected loops, 
a suitable magnetic field should place some inner loops 
in the destructive regime but leave the outer loops super- 
conducting, leading to spatial variation of the supercon- 
ducting energy gap and possible quantized states. 

The question of the existence of quantized states in 
ultrasmall superconductors can be addressed experimen- 
tally. Recent advances in nanofabrication techniques 
have enabled the fabrication of devices with a size needed 
to address this issue. Because the longest ^(0) in thin 
films tends to be on order lOOnm, we fabricated our de- 
vices using ebeam lithography with a single layer PMMA 
resist to achieve smaller feature sizes. All writes were per- 
formed using a Leica EBPG5-HR ebeam writer. Because 
we were interested in obtaining long coherence length 
films, we thermally deposited 20nm of 99.9999% (6N) pu- 
rity aluminum, and defined the devices by a liftoff pro- 
cedure. Devices were fabricated in a quasi- four point 
geometry (see Fig. la for a schematic) where the trans- 
port included a portion of thin Al leads as well as the 
confined superconducting Al square. It should be noted 
that we rely on the difference in the widths of the leads 
and square, rather than normal leads, to provide the con- 
fining geometry. After fabrication we imaged each device 
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FIG. 1: (Color Online) a) Schematic of Al devices used in 
this study. The width and the length of the current and volt- 
age leads to the square vary between devices and connect to 
large Al contact pads; b) Resistance (R) vs. temperature (T), 
R{T), at zero magnetic field (H = 0) showing a superconduct- 
ing transition temperature (Tc) of 1.20 K. Inset: Scanning 
electron microscope (SEM) image of this ^-^500 nm square 
sample. Note that the wavy contact lines were due to an old 
filament installed on the Leica EBPG5-HR e-beam writer and 
the high doses used for writing this device. The organic con- 
taminants shown at the bottom center of the image should 
not affect the device; c) Resistance vs. magnetic field applied 
perpendicular to the substrate surface, R{H). The high and 
low critical fields correspond to those of the thin leads and 
the square, respectively; d) Phase diagram constructed from 
R{H) and R{T) data. Arrows indicate the magnetic fields 
predicted for the transitions of few vortex states in a circular 
disk of the same area as the square (see text) . 



using a scanning electron microscope (SEM) to determine 
the geometry. We constructed critical field verses tem- 
perature phase diagrams of each sample from resistivity 
measurements taken via a dc technique in an RF filtered 
dilution refrigerator. 

The phase diagram in H-T space can be obtained by 
searching for the appropriate magnetic field/temperature 
values for which the sample resistance was kept at a fixed 
value, while the temperature and the magnetic field were 
swept. Because the sample has multiple transitions, due 
to the differing sizes of the leads and square, measure- 
ments of resistance vs. magnetic field at selected tem- 
peratures were used to find the resistance value corre- 
sponding to the superconducting transition of the square. 
This allowed us to construct the H-T phase diagram for 
each of the superconducting squares from our transport 
measurements. Note the convenience of this technique 
in that it allows the determination of the phase diagram 
even when one or two leads are lost by taking into ac- 



count the fixed resistance offset of the leads. In Figure 
1, we show results for Sample A that consists of a square 
roughly 534 nm across and two 94 nm wide leads. In ad- 
dition, two voltage leads are located 600 nm away from 
the square (not shown). The sample has a superconduct- 
ing transition temperature (Tc) of 1.20 K in zero field, 
close to the bulk Tc for Al. To estimate the values of ^(0) 
and the penetration depth A(0), we measured both the 
perpendicular critical field Hcx, using a bulk film from 
the same run, and the parallel critical field Hc|| using 
the thin leads connecting to the square, estimating A(0) 
using [15] 
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where $o = ff is the flux quantum, d is the apparent 
thickness of the parallel portion and He is the thermo- 
dynamic critical field. Combining these two equations 
we can estimate ^(T) by simply measuring Hc|| valid for 
both type I and type II superconductors 
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If our films are type II superconductors, Hcx= Hc2 j^ He 
implying that we cannot accurately estimate A however 
we can estimate ^(T) from the measured Hc2 by 
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We can calculate both A(0) and ^(0) by taking into ac- 
count the temperature dependance ^(T)^ = ^(0)^— j^ 

and A(T)^ = A(0)^ tJ't_ - From the experimental value of 

Tc 

Hc± = 175 G and Hc|| = 1166 G at 100 mK, we obtained 
a ^(0) of 112 nm from the H^n and estimate A(0) of 128 
nm, noting that this implies the sample is type II and 
the value of A(0) from this estimate is not valid. Since 
the film is type II we also estimated ^(0) from Hc± to be 
131nm. It should be pointed out that in the above esti- 
mate was taken for a film evaporated at the same time 
as the squares, so should be similar to ^(0) for the sam- 
ples, while the Hc|| estimation for ^(0) was taken for the 
thin leads connected to the square which should better 
represent the properties of the film at the square. Nev- 
ertheless, because the values of ^(0) estimated here are 
used only to characterize the samples, and are not used 
in any calculation, the rough estimates do not have any 
serious consequence on our analysis of the results. For 
the sample shown in Fig. 1, L/^(0) = 4.79, where L is 
the size of the square. In Fig. Id, we show the H-T phase 
diagram constructed from R{H) measurements taken at 
various temperatures. 



Several interesting features are seen in the phase 
boundary shown in Fig. Id indicating possible transi- 
tions between different states. In addition to a promi- 
nent "bump" near 70 G, two relatively subtle "kinks" 
are also seen near 145 and 270 G, respectively. In order 
to understand the physical origin of these features, it is 
useful to examine the expectations from G-L theory for 
the modulation of Tc as a function of applied magnetic 
field. For disks larger than a critical size the solution 
to the G-L equation features eigenvalues of angular mo- 
mentum corresponding to the mesoscopic disk admitting 
vortices. While the analytical solution for squares is not 
available, numerically derived phase diagrams for vortex 
decoration exist in the literature [16], but don't address 
the low field giant vortex transitions. However the single 
vortex and giant vortex solutions are expected to occur 
at similar fields as in a disk of equal area. We estimated 
the transition fields for n-vortex states for a disk with 
area equal to that of the measured square based on the 
analytic solutions obtained previously [17|. For disks of 
this area, the expected n=0 to n=l transition should oc- 
cur at 145 G and the n=l to n=2 transition should occur 
at 260 G. Therefore the two "kinks" found in the phase 
diagram near these magnetic fields must be associated 
with the few- vortex states of the square, similar to pre- 
vious measurements [7 . On the other hand, no feature is 
expected near 70 G within the G-L theory. The physical 
origin of this "bump" is not understood. 

The observation of few-vortex states in moderately 
confined Al squares is reasonable, because the sample 
size is larger than the size of a vortex. However as the 
sample size is reduced to be comparable with the size of 
the normal core of a vortex, roughly twice ^(T), this vor- 
tex state should cease. Indeed, the solution of the G-L 
equation for a superconducting disk indicates that when 
the radius falls below a critical value of 1.319^ [14 , the 
disk remains vortex free for all magnetic fields below the 
intrinsic critical field. It is important to note that this is 
the regime where G-L theory, which can describe physics 
at a length scale larger then ^(0) (all physical quantities 
in the G-L theory are averaged over this length), starts to 
fail. Therefore how the superconducting order parameter 
responds to an applied magnetic field can not be treated 
self-consistently by G-L theory in this regime. Experi- 
mentally, to approach this regime either the sample size 
needs to be reduced, or the coherence length increased, 
or both. We show in Fig. 2, a 546 nm square prepared 
with C(0)= 209 nm, corresponding to L/^(0)=2.61, which 
has no features identifiable with few- vortex states in it's 
phase diagram (Fig. 2a), indicating that it was in the 
n=0 state for the whole temperature range probed. Sim- 
ilarly for a 250 nm square with ^(0)=102 nm, correspond- 
ing to L/^(0)=2.38, no vortex state was observed. These 
observations suggest that the Al square is too small to 
support a single vortex, as expected. 

We reduced the relative size of the Al squares to ^(0) 
to further understand the effect of confinement on the 
Cooper pairs. In Fig. 3, the phase diagrams of two Al 
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FIG. 2: (Color Online) Critical field versus temperature for 
two squares with their sizes and estimated ^(0) indicated. 
Insets show SEM images of each sample. 
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FIG. 3: (Color Online) Critical field versus temperature for 
a two ultrasmall Al squares, a nominally 240 nm square with 
^(0) = 110 nm (a) and a nominally 140 nm square with a ^(0) 
= 116 nm. The phase diagram for Sample D was obtained by 
three-point measurements, and Sample E by two-point mea- 
surements. Insets show SEM images of the samples with the 
relevant length scales. 



squares featuring sizes of 240 nm and 138 nm, a zero field 
Tc of 1.23 K and 1.28 K and a ^(0) of 112 nm and 117 
nm, respectively, are shown. Given that L/^(0)=2.13 and 
L/^(0)=1.18, a smooth phase boundary similar to those 
shown in Fig. 2 is expected from G-L theory. Surpris- 
ingly, sharp steps are seen in the phase diagrams of these 
two samples. In Fig. 3a, a sharp rise in critical field He is 
seen as the temperature is lowered to below 0.9 K and a 
less distinct rise at 1.23 K (near Tc) as well. For the sam- 
ple shown Fig. 3b featuring even stronger confinement 
than that shown in Fig. 3a, two pronounced, very sharp 
steps are seen, occurring at T = 1 K and 0.6 K. The 
systematic behavior found in these two samples suggests 
strongly that these step features are intrinsic features of 
the Al devices. 

The presence of measurement leads in our sample, how- 
ever, makes even the smallest filled square we studied 
an open system. As a result, we treat the square as 
well as the leads as a whole system in our attempt to 
obtain insight into the origin of the steps seen in the 
phase diagram. For Sample E, ^(0) of the superconduct- 
ing leads is 116 nm and the square size is 138 nm, making 
it reasonable to treat both the square and the leads as 
a quasi one-dimentional (ID) wire featuring two widths 
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FIG. 4: (Color Online) ^(T)(l-T/Tc)^/^ plotted as a function 
of temperature for Sample E with Tc = 1.28K. Inset shows 
large scale SEM image of this device, note the different width 
leads. Apparent streaking in the image due to charging of the 
substrate. 



and branches associated with voltage leads. In the zeroth 
order approximation, we can ignore the complications of 
these leads, and make use of Eq. 3 to relate measured 
Hc|| of the square to ^(T) within the G-L theory even 
though the exact value of ^(0) so deduced may not be 
accurate. Noting that ^(T) - (1-T/Tc)"^/^ we can plot 
^(T)(l - T/TcY^'^ vs. T, as shown in Fig. 4 for sample 
E, showing two strikingly fiat steps. Note that while the 
functional form ^(0)=^(T)(1 - T/Tc)^^^ is only strictly 
valid near Tc the flat steps far from Tc still signify a 
change in ^(0). The implications of Fig. 4 is that the 
quantized states seen in Fig. 3 may be associated with 
different values of ^(0). Essentially, as the temperature 
is lowered to below Tc, the sample starts off with a very 
long intrinsic coherence length that in fact diverges at Tc, 
which reduces to a ^(0) ~ 800 nm near 1 K, then jumps 
to a stable state with ^(0) ~ 408 nm and again at 0.55 K 
to another state with ^(0) ~ 150 nm. Interestingly, the 
same plot of the data for Sample D was found to result in 
very similar results, even though the size of the square is 
much larger than ^(0), making the above argument not 
applicable. 

This suggests that the observed steps in the phase dia- 
gram might be associated with geometrical constraint ef- 
fects of superconductivity in ultrasmall Al squares. Con- 
sider now a quasi ID sample that has exactly the same 
geometry as the experimental sample except that the 
square width is reduced to be the same as the measure- 
ment leads. Obviously no geometrical confinement effects 
are expected. In this limit, in zero field, the supercon- 
ducting order parameter will remain constant throughout 
the sample, with the exception that it may be enhanced 



near the nodes of voltage leads [T8^. With the addition 
of a moderate magnetic field the order parameter should 
remain constant except over a length scale of ^(T) near 
the bulk leads where it will be smoothly varying to ac- 
commodate them being driven normal. As the tempera- 
ture increases, ^(T) increases and diverges at Tc but no 
abrupt changes are expected. Now if we start to enlarge 
a section of the wire gradually and eventually recover the 
original sample geometry featuring a filled square, the or- 
der parameter is expected to evolve correspondingly, es- 
pecially once the enlarged section becomes thick enough 
to provide a "trap" for Cooper pairs. The square and the 
rest of quasi ID system should still feature a well-defined 
order parameter smoothly varying near the joint between 
the quasi ID part of the wire and the square. Since the 
square can not host a vortex, the observed quantization 
could be due to the quantum-size effects of Copper pairs. 
With ^(0) being essentially the size of the Cooper pair, 
quantized values of ^(0) could be associated with quanti- 
zation of the Cooper pair size. Our experimental results 
may suggest that the trapped superconducting conden- 
sate adopts different patterns of the superconducting or- 
der parameter characterized by different values of ^(0) as 
the temperature is raised. 

This analysis, that Cooper pairs can exhibit different 
patterns having different values of ^(0), assumes that our 
geometry can "trap" Cooper pairs, and that these fea- 
tures seen in the data don't arise from some interplay 
between the different width leads in our devices. Our ge- 
ometry could even be insufficient to constrain the Cooper 
pairs. In the smallest samples, the ratio between the 
width of the thin leads, and the size of the square is only 
~30% meaning that the squares cannot fully confine the 
Cooper pairs, and they will leak into the leads. This 
means our observations could be due to a consequence of 
our experimental design rather than effects originating 
from confining Cooper pairs. In addition, our experi- 
mental design complicates the analysis because the two 
thin leads are expected to possess enhanced condensation 
energy, and thus the critical field as well[12 . Near Tc the 
coherence length diverges, and thus the sample also be- 
comes coupled to the wider leads which have a reduced 
He decreasing the observed He. Noting that samples C, 
D and E were not four-point resistivity measurements, 
therefore the phase diagram constructed at a particu- 
lar resistance value may have included the resistance of 
the sample as well as the resistance of the leads. With 
all this in mind, we propose to explain the nature of 
the observed features using proximity effect between the 
different width lead sections having different He along 
the device, see inset of Fig. 4. However, the absence 
of the "quantized" behavior in sample C which differs 
from sample D only in ^(0) indicates that our observed 
features are not explained by a simple interplay between 
the leads and square. More experiments, with better con- 
trolled geometries and possibly tunneling electrodes are 
needed to resolve these issues. 
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